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Abstract 

We consider the asymptotics of flat, radiation-dominated isotropic universes in 
four-dimensional theories with quadratic curvature corrections which may arise 
when contributions related to the string parameter a' are switched on. We show 
that all such universes are singular initially, and calculate all early time asymp- 
totics using the method of asymptotic splittings. We also examine the late time 
asymptotic behaviour of these models and show that there are no solutions which 
diverge as e* 2 , the known radiation solution of general relativity is essentially the 
only late time asymptotic possibility in these models. 



1 Introduction 

A marked feature of theoretical cosmology during the last few decades has been the 
proliferation of different dynamical theories of gravity that are used to tackle the basic 
cosmological issues. With the advent of string theory and its subsequent generalizations, 
these theories acquired a status similar to the shafts and knobs of J. L. Synge's allegor- 
ical box, an unopened (perhaps unopenable) box in which a vast number of shafts are 
connected to one another by many complicated laws, making it objectively impossible to 
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explain what really is inside. There is, it is true, an appealing commonsense in this new 
picture and, although the promise of hard rationality is not yet fulfilled, one hopes that 
in time some new forms will develop and the whole picture will come into better focus. 

Another feature common to many string theory effective actions is the presence of 
higher order curvature corrections, the importance of which cannot be neglected unless 
the basic string parameter a' is small. The generic form of the higher order curvature 
terms present in these effective string actions is that of a typical Gauss-Bonnet combi- 
nation, cf. [U[2]- Such theories admit constant curvature vacua even in the absence of a 
cosmological constant, many of which are obviously stable [3]. 

The pertinent issue of the structure and nature of cosmological singularities becomes 
of interest in this new context for at least two reasons. First, it is important to know 
how different asymptotic behaviours can be compared to the known situation in general 
relativity, for instance whether or not simple isotropic solutions act as past attractors to 
more general homogeneous universes in quadratic theory, cf. [H El [6] . Secondly, in the 
context of string theory the basic problem becomes that of understanding the structure 
of cosmological singularities and how the latter could be resolved using the notion of 
string dualities, cf. [7j. 

In this note we face the more humble task of analyzing the asymptotic form of certain 
simple cosmological solutions to the general quadratic theory in four dimensions. Our 
analysis includes both late and early time asymptotics. We shall show that contrary 
to claims usually made in the literature of such models (such claims were first made as 
early as 1969, cf. [8]), one cannot conclude that solutions (which are regular initially) 
fail to approach the general relativistic solution at late times. In fact, we show that there 
can be no regular solutions at early times and by applying the method of asymptotic 
splittings of [9] , we further prove that there is an open set of initial conditions for which 
the general solution blows up at the big bang (collapse) singularity. 

The plan of this short paper is as follows. In the next Section we write down the basic 
equations of the quadratic theory which we consider further in later Sections. We then 
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show that their late time asymptotics cannot be of the divergent form e* , contrary to 
known claims found in the literature. In Section 3, we apply the method of asymptotic 
splittings to find the form of the general solution of these models near the collapse 
singularity in the past. We may therefore conclude that simple bouncing, flat, radiation 
universes with quadratic curvature corrections do not exist. 

2 Field equations and late time asymptotics 

We start with the action 



<S = / CtatgxdfjLg, d/jg = y/^gdn, (1) 

where £ to tai is the lagrangian density of the general quadratic gravity theory given in the 
fornj] Aotai = C(R) + £ ma tter, with 

C(R) = R + BR 2 + CRic 2 + DRiem 2 , (2) 

where B, C, D are constants. Since in four dimensions we have the Gauss-Bonnet identity, 

5 / R 2 GB dfi g = 0, R 2 GB = R 2 - 4Ric 2 + Riem 2 , (3) 
Jm 4 

in the derivation of the field equations through variation of the action associated with 
(T2]), only terms up to Ric 2 will matter. Therefore the variational derivative of the action 
leads to the following field equations: 



N " T ij = R ij --g ij R + B 



c 4 2 
+ C 



2RW - X -R 2 g^ - 2(g ik gi m - g^g km )V k V m R 



2R ik R\ - \g i] R km R km + V k V k R l] + g l] V k V m R mk - 2V k V J R ki 



(4) 



Below we focus exclusively in spatially flat universes of the form 

ds 2 = dt 2 -b(t) 2 (dx 2 + dy 2 + dz 2 ), (5) 



1 the conventions for the metric and the Ricmann tensor are those of 1 10 



which are radiation dominated (P = p/3). For such spacetimes we have a second useful 
identity, 
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/ (R 2 - 3Ric 2 )dfi g = , (6) 
</m 4 



' M 4 
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which further enables us to include the contribution of the Ric term into the coefficient 
of R 2 , altering only the arbitrary constants. In this case the field equations (jl]) simplify 
as follows: 



2Rffi - l -R 2 g^ - 2(g ik gi m - g^g km )V k V m R 



^ G T\ (7) 



c 4 



where k = QB + 2C . This naturally splits into 00— and %%— components, but only the 
00— component of (|7|) will be used below. This reads: 

b 2 



b 2 K 



2 bb, 2 b&_b^_ b^ 
b 2 b 3 b 2 ¥ 



where b\ is a constant defined by 



8nGp b 2 



(from WP° = 0). (9) 



3c 4 6 4 ' 

Note that the Friedmann solution \ / 2bit of general relativity satisfies the above equation. 

In the rest of this Section we shall consider the problem of the late time asymptotics 
of solutions to Eq. (jSJ). We say that a solution b(t) is asymptotic to another solution 
a(t) provided that the following two conditions hold (the first is subdivided): 

(i) Either (1) a(t) is an exact solution of the system, or (2) a(t) is a solution of the 
system (substitution gives = 0) as t — ► oo, 

(ii) b(t) = a(t)[l + g(t)], g(t) — > 0, as t tends to infinity. 

If either of these two conditions is not satisfied, then b(t) cannot be asymptotic to a(t). 

Let us assume, following [8], that Eq. (JSJ) has a solution with a regular minimum at 
the arbitrary time to, (bo = b(t = to) = and bo = b(t = to) ^ 0). We can then expand 
this solution as a Taylor series 

Kt) =bo+ b i (t -to) 2 + (t -to) 3 + • • • , (10) 
2 o 



and substitute this form back to Eq. (jHJ), to see that this restricts the value of the 
constant k to k = (bi/(b b )) 2 > 0. 

Reduction of the order of Eq. (jHJ) can be achieved if we set / = (bb) 3 ^ 2 and £ = 
12~ 3 / 4 6 3 , so as to obtain the following second order differential equation: 

r-^(r 1/3 -^r 5/3 ) = o, (n) 

where the derivative is with respect to £. For large values of £, an exact asymptotic 
solution of (fTTT) . in the sense of condition (i2) above, is given by 

/ 4 \ 3/4 

/~ — £(ln0 3/4 , for large £, (12) 



3k J 

which in terms of the variables b and t has the form 

b(t) ~ i2V4 e (*-*d) a /i2". (13) 

However, we stress that this is no more an asymptotic solution of the original equation 
(JHJ) (as it is, for instance, claimed in [H]). In fact, a simple substitution of ffTBl into the 
left hand side of Eq. (jHJ) gives the result: 

1 36k' 

which never goes over to zero as t tends to infinity, so that condition (i2) can never 
be satisfied. In other words, an asymptotic solution of Eq. ( fTTT) does not necessarily 
translate into an asymptotic solution of Eq. (jHJ), and this is the case for f[T2"J) . These 
results lead us to conclude that there exist no solutions of the field equation (JTTJ) which 
have late time asymptotics of the form (JTBl . 

One candidate solution (leading to a well-posed asymptotic problem for late or early 
times) that satisfies condition (i) of the definition and so qualifies to perturb is the 
radiation solution proportional to \/t as we already noted after Eq. Qj. The problem 
of the late (as well as early) time asymptotics of solutions to the field equations of the 
quadratic theories was taken up in Refs. [U [S] using a detailed perturbation analysis 
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of the FRW radiation solutions (flat and non-flat) of the form a(t) — (t — at 2 ) 1 / 2 (in 
our present work we focus in the flat case, a = 0). The conclusion is that, although 
the non-flat radiation solutions are generically unstable, in the case of a flat universe 
there is a parameter region in which all late time solutions are stable with respect to 
perturbations, the latter decaying as g (t) ~ t^+r 3 / 4 , (cf., e.g., [5], Eq. (21)). Therefore 
in such models all relevant solutions of our higher order gravity theory asymptotically 
approach the flat, radiation solution of general relativity. In view of the non-existence 
of an e* 2 late time asymptotic shown currently, these perturbation conclusions suggest 
that the flat, radiation form \ft is a unique late time asymptotic solution in the category 
considered here. 



3 Early asymptotic splittings 

We now move on to perform a local asymptotic analysis in order to find the general 
behaviour of the solutions of Eq. (jSJ) near the initial singularity. This analysis is based 
on the use of the method of asymptotic splittings expounded in Ref . [9] , we follow their 
notation closely. As a first step, after setting b = x, b = y and b = z, Eq. (JHI) can be 
written as a dynamical system of the form x = f (x), x = (x, y, z): 

V b i yz ^z 2 3y 3 

x = y, y = z, z = — -- + + 

2k 2nyx z x 2y 2x z 

If a = (a, (3, 7), and p = (p, q, r), we denote by x(r) the solution 

x(r) = ar p = (ar p , {3t\ 7 r r ), (15) 

and by direct substitution in our system (I14p . we look for the possible scale invariant 
solutions of this for There are two possible combinations. The first, which is the 
most interesting one, has dominant part 

f (0) =^,--^ + ^), (16) 
V 2y x 2xV V ; 



2 The vector field f is called scale invariant if f (ar p ) = r p 1 f (a) 
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while the subdominant part reads 

2/q/x 2 2 k 



• sub 



°'°>-^b + 7r:), (17) 



with f = f^ ) + f sub . The dominant balance (of order 3) turns out to be 

<^N'-i)44-!)). < is > 

where a is an arbitrary constant. We recognize this as being what we want, but we 
are not done yet. We now calculate the eigenvalues of an important matrix, called the 
K-matrix, which signify the places in a series expansion of the solution around the finite 
time singularity where arbitrary constants appear, offering thus a clue as to how general 
the found solution is. The Kowalevskaya exponents for this particular decomposition, 
eigenvalues of the matrix /C = Df(a) — diag(p), are { — 1,0,3/2} with corresponding 
eigenvectors {(4, —2, 3), (4, 2, —1), (1, 2, 2)}. The arbitrariness coming from the coeffi- 
cient a in the dominant balance reflects the fact that one of the dominant exponents is 
zero with multiplicity one. 

Keeping with the method of asymptotic splittings [9], we proceed to construct series 
expansions which are local solutions around movable singularities. In our particular 
problem, the expansion around the singularity turns out to be a Puiseux series of the 
form 



z(*) = X) c "(*-*°) i+5 ' y(t) = ^c 2i (t-t )^, z(t) = Y,C3i(t-t )^ } (19) 

i=0 i=0 i=0 

where t is arbitrary and c w = a, c 2 o = a/2, c 30 = —a/4. For these series expansions to 
be valid the compatibility condition 

( -2ci 3 + c 23 ^ 
(1,2,2)- - C2 3+C33 = 0, (20) 

\ ~\ c Vi + | c 23 — c 33 / 

must be satisfied. Substitution of Eq. ffl9|) into Eq. (fT4"|) leads to recursion relations 
that determine the unknowns Cu, Cgj, term by term. After verifying that Eq. (|2"0"|) is 



indeed true, we may write the final series expansion corresponding to the balance ( ITSlh 
It is: 

x{t) =a(t- t )2 + c 13 (t - t ) 2 + OA o 1 (t - to) 1 + • • • • (21) 

The series expansions for y(t) and z(t) are given by the first and second time derivatives 
of the above expressions respectively. 

Our series (ED) has three arbitrary constants, a,ci3,to (the last corresponding to the 
arbitrary position of the singularity) and is therefore a local expansion of the general 
solution around the movable singularity to- Since the leading order coefficients can be 
taken to be real, by a theorem of Goriely-Hyde [TT], we conclude that there is an open 
set of initial conditions for which the general solution blows up at the finite time (initial) 
singularity at t . Finally, we observe that near the initial singularity, all flat, radiation 
solutions of the quadratic gravity theory considered here are Friedmann-like regardless 
of the sign of the R 2 coefficient, while away from the singularity they strongly diverge 
from such forms. This proves the stability of our solutions in the neighborhood of the 
singularity. 

The second possibility of a scale invariant solution has dominant part 

\ ' ' 2nyx 2 2y x 2x 2 ^) ' ^ ^ 



f sub = fo,0,^-) . (23) 



and subdominant part 

f sub = fo.o. 

2k 

The dominant balance (of order 2) at the singularity turns out to be 

{(a, p)} = {((«, a, 0), (1, 0, -1)), ((-«, -a, 0), (1, 0, -1))}, (24) 

where a a constant defined by a = (b\/'$K) 1 l A E K.. The Kowalevskaya exponents are 
{ — 1, VQ, — v^}, therefore this decomposition leads to a particular solution (two arbitrary 
constants) whose dominant term near the initial singularity is b(t) ~ (t— to) and therefore 
does not provide any bounce solutions. 
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4 Conclusions 



We have considered the behaviour of flat, radiation-dominated solutions to the general 
quadratic theory of gravity in four dimensions. This theory is supposed to arise as part 
of an effective string theory action truncated to first order in a'. We found that in these 
models there are no bounce solutions, all have a collapse, initial, isotropic singularity 
and we have given explicit forms of the approach to the initial finite time singularity. 
However, our results suggest that it is not possible to inhibit the existence of late time 
behaviour similar to the observed, even though any such solution must necessarily be 
singular initially. 

There are many ways to extend these results. An obvious one is to take into account 
the dilaton dependence present in the full string action (cf. pQ, Eq. (4.3)), in other words 
to consider the present problem in the string frame (and not in the Einstein conformal 
frame of the present paper). We also wish to see how our results are altered (or not 
altered) when we pass on to models with curvature, especially near the singularity. 
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